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Abstract. This article deals with the numerical approximation of effective coefficients 
in stochastic homogenization of discrete linear elliptic equations. The originality of 
this work is the use of a well-known abstract spectral representation formula to design 
and analyze effective and computable approximations of the homogenized coefficients. 
In particular, we show that information on the edge of the spectrum of the generator 
of the environment viewed by the particle projected on the local drift yields bounds 
on the approximation error, and conversely. Combined with results by Otto and the 
first author in low dimension, and results by the second author in high dimension, 
this allows us to prove that for any dimension d > 2, there exists an explicit numerical 
strategy to approximate homogenized coefficients which converges at the rate of the 
central limit theorem. 
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1. Introduction 

We consider a discrete elliptic operator —V* • AV, where V*- and V are the discrete back- 
ward divergence and forward gradient, respectively. For all z S Z, d , A(z) is the diagonal 
matrix whose entries are the conductances uj z z + ei of the edges (z,z + ej) starting at z, 
where {e^i^ii,...^} denotes the canonical basis of The values of the conductances are 
random and their realizations are assumed to be independent and identically distributed. 

Provided that the conductances lie in a compact set of R+, standard homogenization 
results (see for instance [7]) ensure that there exists some deterministic matrix j4hom 
such that the solution operator of the deterministic continuous differential operator —V • 
^homV describes the large scale behavior of the solution operator of the random discrete 
differential operator —V* • j4V almost surely. As a by-product of this homogenization 
result, one obtains a characterization of the homogenized coefficients ^4hom : it is shown 
that for every direction £ £ M. d , there exists a unique scalar field <f> such that V</> is 
stationary, (V0) = (vanishing expectation), which solves the corrector equation 

- V* • A{£ + V(/)) = inZ d , (1.1) 

and normalized by (p(0) = 0. With this corrector, the homogenized coefficients j4hom can 
be characterized as 

£ • AcmC = ((£ + V$ • A(£ + Vfi) . (1.2) 
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From the practical point of view, (|1.2j) is not of immediate interest since the corrector 
equation (jl.ip has to be solved 

• for every realization of the coefficients u, 

• on the whole Z, d . 

Ergodicity allows one to replace the expectation by a spatial average (on increasing do- 
mains) almost surely. To approximate <fi, one usually uses (j)R, the unique solution to 
equation (jl.ip on some large but finite domain Qr = (— R/2, R/2) d , completed by say 
periodic or homogeneous Dirichlet boundary conditions. Yet, the comparison of V^r to 
Vcj) is not obvious since Vc^r and V</> are not "jointly stationary". In order to avoid this 
difficulty, Otto and the first author have used a somewhat different strategy. We have pro- 
ceeded in two steps: we first replace (ft by its standard regularization cft^, unique stationary 
solution to the modified corrector equation 

- V* • A(£ + V</v) = in 7L d 

for some small [i > 0. Then, ^ is replaced by fi^^R, the unique weak solution to 

M<Am?-V* -A^ + V^r) = mQ R nZ d , 
<t>^ R = onZ d \Q R . 

The advantages are twofold: 

• V</> and V<^ are jointly stationary, which is of great help for the analysis, 

• (j)n is accurately approximated by 0„ r on domains of the form Ql = (—L/2,L/2) d 
provided that {R—L)^fjl S> I, due to the exponential decay of the Green's function 
associated with fj, — V* • AV in Z d (see [2]), so that we only focus on <ft^ and not 
4>fj, t R from now on. 

In particular, we may approximate Ah om by the following average 

Z-A^lZ-- I (£ + V</vM(£ + V^)xz(z)dx, 

JQl 

where \L is a smooth mask supported on and of mass one. In [3], we have proved that 
the L 2 -norm of the error in probability takes the form 

( (£ • A^U - £ • Aom£) 2 ) = var [£ • A^tf] + (£ • (A^ - A hom )tf, (1.3) 

where 

s-A^it ■= ((e + v^)-A(e + v^)). 

The first term of the r. h. s. of (|1.3p is stochastic in nature and corresponds to the 
variance of the approximation of the homogenized coefficients, whereas the second term 
is a systematic deterministic error related to the fact that we have modified the corrector 
equation. 

In [3J, we have proved that the stochastic error depends on the dimension and has the 
scaling of the central limit theorem (in other words the energy density of the corrector 
behaves as if it were independent from site to site): there exists q depending only on the 
ellipticity constants a, j3 such that 

L~ x hi q u if d = 2, 

L-V* if d>2. (L4) 
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The systematic error has been identified in [3]. It also depends on the dimension for d < 5, 
but saturates at d = 5: there exists q depending only on the ellipticity constants a, f3 such 
that 

/xlnV^ 1 if d = 2, 
/z 3 / 2 if d = 3, 



^4/^,1 ^4homl 5n 



/i 2 In /J" 1 if d = 4, (L5) 
/i 2 if d>4. 

These two estimates are optimal (up to some possible logarithmic corrections for d = 2). 
In order to use (fr^ R as a proxy for 0„ on Q^, at first order we may take ~ L 2 ~ i? 2 . 
Hence, the stochastic error dominates up to d = 8, so that the convergence rate of the 
numerical strategy is optimal (it coincides with the central limit theorem scaling, which 
is an upper bound): 

L _1 ln«Z if d = 2, 
L -d/2 j£ 2 <d<8. 



^■A L -, XL i-i-A hora if) 1/2 < 



Yet, for d > 8, the systematic error dominates and the numerical strategy is not optimal 
any longer: 

(^i-M,^-^hom£) 2 ) 1/2 < £~ 4 if 8<d. 

The aim of this paper is to introduce new formulas for the approximation of A^ ora using 
the modified corrector (possibly with different /x's) in order to reduce the systematic 
error. In early and seminal papers on stochastic homogenization (for instance |9j and [5]), 
spectral analysis has been used to prove uniqueness of correctors, and devise a spectral 
representation formula for A^ ora . In particular, denoting by — C the generator of the 
environment viewed by the particle, and by eg its spectral measure projected on the local 
drift d = V* • A£ (see Section [3]) , we have 



As noticed by the second author in [8], A^i can also be written in terms of the spectral 
measure e (see Section [2] for details): 

Z-A^d = (Z-AO- [ ^±^de (A) 

= i-A hom ^ + ^x 2 I — — — 2de (A). 

Jr+ + xy 

The key idea of the present paper is to use this spectral representation in order to design 
approximations of Ah om at an abstract level first, and then go back to physical space and 
obtain formulas in terms of the modified correctors (p^. We shall actually introduce, for 
every integer k > 1, an approximation A^ of Ai lom defined in terms of . . . , ^-u, 
and prove that, up to logarithmic corrections, the difference |Ah om — ^u,fe| is bounded in 
our discrete stochastic setting by 

min{2fc,d/2} -r j < n 
^min{2fc,max(3,d/2-3)} jf d > 6 

(see Theorem [3] for a more precise statement). The systematic error associated with the 
new approximations can be made of a higher order than (jl.5j) as soon as d > 4. The proof 
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of these estimates relies on the observation that the systematic error is controlled by the 
edge of the spectrum e ((0, //)). In turn, the systematic error also controls the edge of the 
spectrum (see Theorem |4] for a precise statement), so that estimating the systematic error 
is equivalent to quantifying e ((0,/x)). 

As we shall also prove, the variance estimate (|1.4|) is unchanged if A„i is replaced by 
A„k for all k > 1. In particular, if we keep fi^ 1 ~ L 2 , we obtain a numerical strategy 
whose convergence rate is optimal with respect to the central limit theorem scaling in 
the stochastic case, for any d > 2. This improves and completes for d > 8 the series 
of papers [3j HI [2] by Otto and the first author on quantitative estimates in stochastic 
homogenization of discrete elliptic equations. In turn, we also obtain "optimal" bounds 
on e ((0, (J,)) up to d = 6 (see Theorem [5]), thus improving the corresponding results of the 
second author in [8]. 

Note however that the bounds (jl.6p are not yet optimal: the systematic error is expected to 
behave as ^{2k,d/2} in any 

dimension (up to logarithmic corrections), see also Remark [1] 
for the equivalent statement in terms of the edge of the spectrum. We wish to address 
this issue in a future work. 

The article is organized as follows. Although the main focus of this work is on stochastic 
homogenization of discrete elliptic equations, we first describe the strategy on the elemen- 
tary case of periodic homogenization of continuous elliptic equations in Section [2] (this new 
strategy may indeed be valuable to numerical homogenization methods, see in particular 
PQ for related issues). We introduce the spectral decomposition formula for the homog- 
enized coefficients. The binomial formula then provides with natural approximations of 
the homogenized coefficients in terms of the associated spectral measure. We conclude 
the section by rewriting these formulas in physical space using solutions to the modi- 
fied corrector equation, which yields new computable approximations of the homogenized 
coefficients. In particular, this generalizes the method introduced in [T] and makes the 
systematic error decay arbitrarily fast. Some numerical tests displayed in Appendix [B] 
illustrate the sharpness of the analysis. 

We turn to the core of this article in Section O the stochastic homogenization of discrete 
elliptic equations. We first recall the spectral decomposition of the generator of the en- 
vironment viewed by the particle. The algebra is the same as in the continuous periodic 
case, so that the formulas we obtain in Section [2] adapt mutatis mutandis to the discrete 
stochastic case. Yet, the error analysis is more subtle. We show that the asymptotic 
behavior of the systematic error is driven by the behavior of the edge of the spectrum of 
the generator. Using results of [8] in high dimension, and results in the spirit of [4] (see 
Lemma [5] and Appendix [AJ in low dimension, we obtain estimates on the edge of this 
spectrum, which show that the systematic error is effectively reduced in high dimensions 
(although our bounds are not optimal when d > 6). We then note that the variance esti- 
mates derived in [3] also hold for these approximations, thus concluding the error analysis 
of the numerical strategy. 

We will make use of the following notation: 

• d > 2 is the dimension; 

• In the discrete case, f^ d dx denotes the sum over x € and f D dx denotes the 
sum over x G Z d such that x € D, D open subset of M. d ; 
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• (•) is the average in the periodic case, and the expectation in the stochastic case; 

• var [•] is the variance in the stochastic case; 

• < and > stand for < and > up to a multiplicative constant which only depends 
on the dimension d and the constants a, (3 (the ellipticity constants of the matrix 
A, see Definitions Q] and [5]) if not otherwise stated; 

• when both < and > hold, we simply write ~; 

• we use 3> instead of > when the multiplicative constant is (much) larger than 1; 

• (ei, . . . , e^) denotes the canonical basis of M. d . 

2. The continuous periodic case 

Definition 1. Let A : W 1 — > A^rf(M) be a Q-periodic symmetric diffusion matrix which is 
uniformly continuous and coercive with constants f3 > a > 0: for almost all x G Q and 
all £ G M. d , \A£\ < and £ • A£ > a|£| 2 . The associated homogenized matrix ^4hom is 
characterized for all £ G M. d by 

t-A hom £ = ((^ + v^)-A(e + V0)) 

= / (£ + V4>) ■ A(£ + V<j>)dx, 

JQ 

where (•) denotes the average on the periodic cell Q, and (f> is the unique Q-periodic weak 
solution to 

- V-A{i + V<j>) = (2.1) 

with zero average (</>) = 0. 

Let us define £ : H^ cr (Q) x H^ er (Q) — > R + , (ip,x) ^ Jq ' AVxdx the bilinear form 
associated with A. We call the quadratic form ip — > £(ip, ip) the Dirichlet form. One may 
write the homogenized matrix as 

e- Aom£ = -AO -£((/), 4>). (2.2) 
Indeed, the weak formulation of (|2.ip implies 

/ • A(£ + V0)dx = 0, (2.3) 
JQ 

and therefore 

/ t-AVcpdx = [ V<t>-A£dx W - f Vcp-AVcpdx = -£(</>, <f>). 
JQ JQ JQ 

The objective of this section is to use a spectral decomposition to design approximations 
for £ (0, (f>). 

2.1. Spectral decomposition. 

Definition 2. Let A be as in Definition [TJ We let £ _1 denote the inverse of the el- 
liptic operator C = — V • AV with periodic boundary conditions on Lq(Q) = {v G 
L 2 (Q) | Jqv{x)oIx = 0}. It is a well-defined compact operator by generalized Poincare's 
inequality, Riesz's, and Rellich's theorems. 
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By Hilbert-Schmidt's theorem, there exist an orthonormal basis {tp%}i>o of Lq(Q) and 
positive eigenvalues {Aj}j>o (in increasing order) such that for all i > 0, = j-i^i- 

By definition, ipi G H^ CT (Q). Setting tpo = 1 and Ao = 0, one may then characterize 
H leAQ) as 

#per(<2) = \ U = ^2 X^ 1 + Ai ) a * 2 < °° f • 

I iGN iGN J 

By Riesz's representation theorem, this also implies for the dual H~^.(Q) of H^ CT (Q): 

h-^aq) = |/ = E^ ErTA" <oo j- ( 2 - 4 ) 

l ieN iGN * J 

Hence, for all / G H^JQ) such that (f,l) H -i H1 = 0, the unique weak solution u G 

-V • iVu = / 

is given by 

u= E x 

iGN\{0} 1 



For all / G H^{Q) such that (/, 1) = 0, we define the spectral measure ej of £ 

projected on / by 

e / = E ^H^,Hi eI S ^ ( 2 - 5 ) 

iGN 

where <5;^ is the Dirac mass on Aj. The above characterizations of -ff 1 and then allow 
us to give a mathematical meaning to the formal functional calculus 



</,*(£)(/)> ,^ cr = f *(A)d e/ (A), 

p per JR+ 



for every continuous function ^ : [0, +oo) —> R such that A^A) < 1 as A — )• oo. 

We are now in position to express the Dirichlet form of the corrector (f> in terms of the 
spectral measure projected on the "local drift" 

9 := V • A£ G H-i(Q). (2.6) 

In particular, 

£{(/),(/)) = (C<f>,<t>) H -i H i er 

^de (A). (2.7) 

Let us then turn to the approximation of A^ om used in pQ, that is 

£ • := ((£ + V^) • A(£ + V<^)> , (2.8) 

where ^ G H^ er {Q) is the unique weak solution to the modified corrector equation 

- V • A(£ + V<^) = 0, (2.9) 
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that we more compactly write as (/j, + C)^^ = t). In this case, the weak formulation of 
the equation implies 



[ V^-A(e + V^)dx = -fi [ <j>ldx, (2.10) 
JQ JQ 

so that the defining formula for A„ turns into 



£ ■ Arf = (£ ■ AO - Sfa, ^) - 2fi {<%) . (2.11) 
Proceeding as above, we rewrite the last two terms of the r. h. s. of (|2.11|) as 

deo(A), (2.12) 



(M + A) 



and 



perj-'-'per 

de 9 (A). (2.13) 



The combination of ([22]), (pTTj) . (pTT) . (pJ2|) & (f2TTTfl) allows us to express the difference 
between A^ and A^om in terms of the spectral measure of C projected on the local drift 
D, as observed in [SI Addendum]: 



e • (a, - A hom )e = f(0 l 0)-5(0 M ,^)-2/*<^; 

1 A 2/x 



A (/i + A) 2 Qu + A) 2 
deo(A). 



de a (A) 



'r+ A(/i + A) 2 

Not only does this identity suggest that \A„ — ^4hom| ~ (as proved by a different 
approach in [I]), but it also gives a strategy to construct approximations of Ah om at any 
order. In particular, for all k 6 N, we write 

\2fe 



Jr+ A(/x + A) 2fe 



and set 



(// + A) 2fc - fi 2k 
R + X(fi + A) 2 

'2k\ f u j \ 2k ~ 1 ~ j 



t-A^Z := ({-AO - j , ?W de >^ 



<^o-e(T)/ r+ t^f^>- (2 ' 14) 
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Note that the only operator which has to be inverted to compute A^k is indeed fi + £, 
and not C, as desired. In addition, this definition implies 

£ • (A^k ~ A hom )£ = / — — --TT2fcde (A), 

Jr+ Kfj> + a) lk 

which suggests that the error is now of order /i 2fe . 

In view of formula (j2.14j) . the effective computation of (/i + £)~ k D is needed in practice 
to obtain A^. This is a big handicap for the numerical method since the numerical 
inversion of fi + C has to be iterated k times, which dramatically magnifies the numerical 
error. Fortunately, one may use a sligthly different approximation of j4hom which avoids 
this drawback, as shown in the following subsection. 

2.2. Abstract approximations. Let us first introduce functions Qn,k, which are defined 
as linear combinations of (p^, ... , </>2 fe - 1 A t (and therefore easily computable) and will serve 
as substitutes for (/i + C)~ k D. 

Definition 3. Let A and Q, C, and Z) be as in Definition [U Definition [21 and (|2.6p . 

respectively. For all \i > 0, the sequence of functions £ H^ cr (Q) is defined by its first 
term 

= ^ = fa + Q' 1 *, ci = l, (2.15) 

and by the induction rule 

= c fc ^ _1 (t) Mi fc - t) 2At ,fc), Cfc+i = f^- + 1^ . (2.16) 

Defined this way, the functions satisfy the following fundamental properties: 

Proposition 1. Lei 0^ be as in Definition^ then for all \i > and k > 1, we have 

cV, fc+ i = {fi + Cy^k, (2.17) 

&>n,k+i = (1 + c k )X>2n,k ~ Cfcf^fc- (2-18) 

Proof. Identity (|2.18p is a direct consequence of (J2.16P & (|2.17p , and we only need to prove 
the latter. We proceed by induction. Let us first check that it is indeed true for k = 1. 
By definition of t>^,i, we have 

(^ + £)cV,i =5, (2.19) 

and as a consequence, 

(/x + £)B 2/1 ,i =3-M 2/1 ,i- ( 2 - 2 °) 
Combining (|2.19p and (|2.20p . one obtains : 

(fi + £)(^,i - t) 2At ,i) = /U D 2j u,i> 

from which it follows that J)„ 2 = + C)~ 1 d2 fJ ,,i- Let us now assume that (|2.17p is satisfied 
at level k > 1. Similarly, we have 

Gu + £)9 M ,fc+l = *> 2/i ,fc, 
(M + Q^^k+l = ^Afi,k — fJ> f 2/i,fc+l- 

Using these equalities, together with the definition (j2. 16|) of D^jt+i, we are led to 

(M + ~ 5 2/i,fc+l) = f 2 /i,fe - ®4n,k + M 5 2p,k+l = A 4 + l) 5 2/i,fc+l, 
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and thus, t> M ,/c+2 = (fJ. + C) 1 t>2^fc+i • □ 
In order to be consistent with (|2.17p . we set fl^o = ^ f° r an /-t > 0. 

We are now in position to define a suitable approximation of ^4hom- The idea is to use the 
identity 

1 (/i + A) 2 (2^ + A) 2 ...(2 fc -y + A) 2 
A A(/i + A) 2 (2/i + A) 2 ...(2 fc -V + A) 2 

in (|2.7p . expand, and take advantage of Proposition [T] to efficiently compute terms of the 
form 0^ = (fi + C)~ 1 (2fi + C)^ 1 . . . (2 fe_1 / u + £) _1 £). This gives rise to the following 
(abstract) approximations of A^ om , and systematic errors: 

Theorem 1. Let A and ^4hom be the Q-periodic diffusion matrix and the associated homog- 
enized diffusion matrix of Definition^ For any fixed £ E M. d such that |£| = 1, we denote 
by e a the spectral measure (|2.5p of £ = — V • AV projected on the local drift Z) = V • 
For a// € N, we let P^ : R x R — >• R 6e i/ie polynomial given by 

P fc (/i, A) = A" 1 ((/x + A) 2 (2^ + A) 2 • • • (2 fc ~ V + A) 2 - 2 fc ( fc "V fc ) , (2.21) 

and for all [i > 0, we define the approximation An^ of ^4hom 

i ■ ^ = « • - / k+ ( , + A)2(2> , + ^. A) (2t - V + A)2 ^(A). (2.22, 
T/ien i/te systematic error satisfies 

/ 1 N 2fe 

< e • (^fc - A hom )£ < 2 fe ( fe " 1 ) <|A| 2 ) U±-\ (2.23) 
Proof. Starting point is the identity 



Z-(A,, k -A hom )( = /_ (x- (A + M )2...(A + 2^V) 2) dg ° (A) 

2fc(fc-l)^2fc 



A(A + /i) 2 ---(A + 2 fe -V) 2 



de (A), 



which is a direct consequence of (12. 2p , (12. 7p , and (|2,22p . From this identity, and using Def- 
inition [21 we infer that the systematic error is smaller than and asymptotically equivalent 
to Cfi 2k (as \x tends to 0), where C > is given by 

C := 2*M / de,(X) = 2^ £ ^g^ - (2.24) 

In order to estimate C via (|2.24p . we compare the spectral gap Ai of C to the spectral gap 
\® of — A on Hp CT (Q). By comparison of the two Dirichlet forms, we have 

Ai > aX° v 

The spectrum of the Laplace operator on H^ CT (Q) is explicitly known, and the spectral gap 
given by \® = 4-7T 2 . Hence, recalling that (Z), 1) = and using the characterization 
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dZH) of ffper(Q), one may bound the r. h. s. of (T2T2D by 




as desired. □ 

2.3. New formulas for the approximation of homogenized coefficients. In this 
subsection, we show how to rewrite the approximations A^ tk of A\ lom introduced in The- 
orem [[] in terms of the modified correctors cft^, 02^, . . . , 4>2 k - 1 ^- We proceed by induction. 

Proposition 2. Let c k be as in Definition^ We define the sequence {rafc,i}fc>i,£e{o,...,ifc-i} 
by ai,o = 1 and the induction rules 

a k+lji = c k a kii - 2 1 ~ k c k a k)i _ 1 for i € {1, k — 1}, 
a k +i,k = —2 1 ~ k Ck a k,k-i- 

Within the assumptions and notation of Theorem GJ the approximations A^ <k of j4hom 
satisfy the formula: for all £ € M. d , 

= ((e + v^)-A(e + v^)) 

k-l 

i=0 i=0 j>i 

where the {^lAi^^i are ^ e m °dificd correctors associated with £ through (j2.9j) . and i/ie co- 
efficients {r?fe,j}fe>i,o<j<fe, a™d { zy fc,i,j}fc>2,o<j<fc,o<j<j are defined by the initial value r?i j0 = 
0, and i/ie induction rules 

T]k+i,i = %,i + (2 fc(fc - 1)+i -2 fc2+1 )al +1 , / O rie{0,£;-1}, 

Vk+l,k — ~ z a k+i,ki 

v k+i,i,k = — 3 x 2 k )a k+ i ji a k+ i :k , 

Vk+i,ij = Vk,ij + (2 fc(fc ~ 1) (2 i + 2*) - 2 k2+2 )a k+l!i a k+1J for j € {0, k — 1}. 
./Vote t/iat {^fc,i,j}fc>i,o<j'<fc,o<i<j does not require further initialization. 
Proof. We proceed in four steps. 
Step 1. Proof that for all > 1, 

£ • ^, fe+ i^ = £ • A^Z - 2 k ^ fJ 2k £(D fltk+1 ,^ k+1 ) - 2 fc2+ V fc+1 (t>l k+l ) . (2.26) 

In order to prove (|2.26|) , we first note that the polynomials P k defined in (|2.2ip satisfy the 
identity 

Pk +1 (ji, A) = (2*> + A) 2 P fc (/i, A) + 2 k ( k ~ 1 hfi 2k + 2 fc2+ V 2fc+1 - (2.27) 
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Hence, formula (I2.22p implies 

P k+1 {fj,,X) 
Gu + A) 2 (2/i + A) 2 ---(2V + A) 
G2B lt ... /" (2 fc // + A) 2 P fc (//,A) 



" M) ~ L (^ + A) 2 (2 M + A) 2 ...(2V + A) 2deo(A) 
/• 2 fc ( fc - 1 )A^ 2fc + 2 fc2+ V 2fc+1 
" /r + (^ + A) 2 (2/i + A) 2 ..-(2V + A) 2de5(A) 
a , a c f 2 fc (^ 1 )A/i 2fc + 2 fc2 +V 2fc+1 ^ 
" ^ ^ L (/, + A) 2 (2^ + A) 2 ...(2^ + A) 2de5(A) - 

From (|2.17p in Proposition [Q we infer that (/i + • • • (2 fc /i + £)~ 1 t) = D^k+i, so that 
the above identity turns into 



jjp 0r ,jj.p Cr 

as desired. 

Step 2. Proof that for all k > 1, 

fc-i 

M fc-1 iWfc = J>k,i&V ( 2 - 28 ) 
i=0 

We proceed by induction, and assume that (|2.28j) holds at step k. The induction rule 
(|2.16p then yields at step k + 1 

= ^ k ~ l Ck{^)^k-^2^k) 

= CfcM* -1 ^* " c fc 2 1 - fc (2/i) fc " 1 5 2At ,fe 
fc-i fe-i 

i=0 i=0 

= c k a k ,o4>ii + f ^Cfc(afe,i - 2 1 ~ fc a fe ,i-i)</>2v J - c k 2 1 ~ k a k ^ l 4> 2kfl , 
^ i=i ' 

so that fi k d^ tk+ i = YM=o a k,i4>2 i k Li as desired. It remains to recall that D^i = </>^ to 
conclude the proof of f|2.28j) . 

Note that 02,0 = 1 and 02,1 = — 1. In particular, 02,0 + 02,1 = and the property 

fc-i 

£ OJM = (2.29) 

i=0 

follows by induction, for all k > 2. 

11 
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Step 3. Proof that for all i,j > 1, 

£(&v>^m) = (e-^)-^(<(e + v^ 2 g-^ + v0 2> )) 

+<(e+v^ v )-^(e+v^ M ) 

-M^" 1 (^(^ + ^ M )> - / u2 j - 1 (fonifon + <hfu)) ■ (2-30) 

We can easily see that (|2.30p holds when i = j. From (|2.8|) and (|2.1ip . we have indeed 
that 

sifovfoj = (e-^)-<(c+v0 2 g^(e+v0 2 g>-2 i+ v(^ v ), (2.31) 

as desired. For general i,j £ N, we have, using (|2.10p first for ^> 2 j M and then for 

= (V^ • A(£ + V0 2JM )) - <V0 2V • A£) 

IP -2V {fo^fa^ - <V0 2> • A(£ + V^)) + (V^ • AV(j> 2itl ) 



^ -2V<^^,)-2V(4,) 

+ <£ • ^) - ((£ + V0 2 g • A(£ + V^ M )> . 

We conclude the proof of (|2.3U|) by changing the roles of i and j. 

Step 4. Proof of (f2T25j) . 

In view of (|2.26p . we have to estimate two terms. We begin with the Dirichlet form: 
inserting (|2.28p in the integral yields 

k k 
i=0 j=0 

We then appeal to ([2.30P to turn this identity into 
M 2fc £(^V,fc+i> <V,fc+i) 

■ fc 2 k k 



j=0 i=0 ,J>i 

fe 

i=0 



2 a fc+lj4 C 2 a fc+ ij) ^2V + ((£ + V0 2V ) • A(£ + V0 2> )) 



i=0 v j=0 

12 
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Taking into account (|2.29p . we finally have 

k k k 

= ~ E Y ^ + 2 J 'H +li ia fe+ i, j ((fopten) - Y & ia k+i,i ■ ( 2 - 32 ) 

We now turn to the last term of the r. h. s. of ()2.26|) and appeal to (|2.28D : 

k k k 

V 2k+1 (fljU+l) = Y M a fe+l,i (^2>) + 2 E Y MOfc+l,iOfc+lj (02*m^2»>) • ( 2 - 33 ) 
i=0 i=0 j>i 

We then prove (|2.25p by induction, recalling that 

a m = ((e + v^)-^(e + v^)). 

Let us assume that (|2T25]) holds at step k > 1. Combined with (f!T32"|) & (p33|) . flgjggp 
turns into 

e • a^^^ = e • ^ fc $ - 2*< fc - v fc £(Vk+i,Wi) - 2 fc2+ v 2fc+1 (^ )fc+1 > 

fc-i fc-ife-i 

+M X] ^ (^2*m) + ^ Y Y Vk ^ (^V^W) 

j=0 i=0 j>i 

i=0 j>i i=0 

-2 fc2+ V ^ (<^ v ) - 2 fc2+ V X] X Ofc+l,i<»fc+l J (02>02^) , 

i=0 i=0 ,j>i 

from which we deduce that ()2.25|) holds at step k + 1. □ 
Proposition [2] yields the following formulas for the first four approximations of Ah m : 

£-A M , 2 £ = ((e + V^)-A(e + V^))-3 / u(^>-2^(^) + 5^(^ 2At ), 

55 4 
e-A^ = (^ + V0 M )-A(^ + V0 M ))- ¥ /i<^)-8/x<^)--/i<0l M ) 

41 22 10 

e-A M e = ((^+v^)-A(e+v^))-^<^)-^<^)-^<^) 

8 . 2 . 1325 ^ . 370 t± ± . 184 ., 



441- \ ~ °r/ 63 - "c^' 63 r 441 

82 /- - V 44 /. 1 V 20 

+ {<P2n<P4n) - g^M W2finy) + go 
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2.4. Complete error estimate. In this subsection, we combine the approximation for- 
mulas with the filtering method used in [T]. The filters are defined as follows. 

Definition 4. A function \ '■ [ — 1 , 1] — > K + is said to be a filter of order p > if 

(i) X GC7P([-i,i])nWP +1 - 0O ((-i,i)), 

(h) Li x(x)dx = 1, 

(iii) x (fc) ("l) = x (fc) (l) = for all k € {0, . . . ,p - 1}. 
The associated mask xl ■ L] d — > M + in dimension d > 1 is then defined for all L > 
by 

d 

Xl(x) := L^HxiL^Xi), 
i=i 

where x = (x\ , . . . , Xd) € 

Let now A and Ah om be as in Definition [TJ For all k > 1, // > 0, p > 0, and R > L > 0, 
we define the approximation A^fe^i of ^4hom a s 

^■^W? := (((£ + V^ B ) ■ A(£ + V^ |fl ))) L 
fe-l fe-i fc-i 

+f J '^2 r )k,i(( ( l ) 2* fJ „R))L + M^^^W^^^ftfl})!, (2-34) 
i=0 i=0 j>i 

where the coefficients 77^ and ffe,i,j are as in Proposition [21 the modified correctors 
are the unique weak solutions in Hq(Qr) to 

2V02V,ii - V • A(£ + V^ >fl ) = 0, 
and denotes the average with mask xl- 

(( h ))L ■= / h(x)xL(x)dx. 

The combination of [H Theorem 1] with Theorem Q] and Proposition [2] then yields 

Theorem 2. Lei d > 2, A and ^4hom ^ e as *n Definition^ k > 1, x be a filter of order 
p > 0, and j4^,fc,_R,L fre i/ie approximation (|2.34|) 0/ £/ie homogenization matrix, where 
R 2 > /i" 1 > R, R > L ~ R ~ i? — L. Then, there exists c > depending only on a, /3 and 
d such that we have 

\A ML - A hom \ < + fi 2k + M ~ 1 / 4 exp {-c^Jl{R - L)) . (2.35) 

In order to illustrate Theorem [21 we provide the results of numerical tests in a periodic 
discrete case in Appendix [Bj They confirm the sharpness of the analysis. 

3. The discrete stochastic case 

We start this section by defining the discrete stochastic model we wish to consider. 

14 
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3.1. Notation and preliminaries. We say that x, y in Z are neighbors, and write x ~ y, 
whenever \y — x\ = 1. This relation turns Z rf into a graph, whose set of (non-oriented) 
edges we will denote by B. We now turn to the definition of the associated diffusion 
coefficients, and their statistics. 

Definition 5 (environment). Let O = [a, /3] B . An element to = (cj e ) ee B of Q, is called an 
environment. With any edge e = (x, y) G B, we associate the conductance to x> y := to e (by 



construction to. 



to 



y..r 



). Let v be a probability measure on [a, /3]. We endow with the 
product probability measure P = z^® B . In other words, if to is distributed according to the 
measure P, then (w e ) eg B are independent random variables of law v. We denote by L 2 (0) 
the set of real square integrable functions on for the measure P, and write (•) for the 
expectation associated with P. 

In the framework of Definition [U we can introduce a notion of stationarity. 

Definition 6 (stationarity). For all z G Z rf , we let 9 Z : VL — > be such that for all to G 
and (x,y) G B, (9 Z to) XtV = oj x+Zty+z . This defines an additive action group {9 z } z£Z d on O 
which preserves the measure P, and is ergodic for P. 

We say that a function / : Q x Z d — > R is stationary if and only if for all x, z G Z d and 
P- almost every to G 0, 

f(x + z,u) = f{x,6 z to). 
In particular, with all / G L 2 (J7), one may associate the stationary function (still denoted 
by /) Z d x Q — > M, (x,to) t— > f(9 x to). In what follows we will not distinguish between 
/ G L 2 (0) and its stationary extension on Z d x $7. 

It remains to define the conductivity matrix on Z rf . 

Definition 7 (conductivity matrix). Let £1, P, and {9 z } z£ i i d De as i n Definitions [5] and [H 
The stationary diffusion matrix A : T, d x — > A^d(M) is defined by 

A(x, to) = diag [w X)X+ej , • • • , w XjX+e J • 

For each to G fi, we may consider the discrete elliptic equation whose operator is 

-V* ■ A(-,to)V, 



where V and V* are defined for all u : Z 

u{x + ei) — u(x) 



V«(i) := 



u(x + e^) — u(x) 



by 



V*w(x) := 



u(x) — u(x — ei) 



u(x) 



e<i) 



(3.1) 



and the backward divergence is denoted by V*-, as usual. The standard stochastic ho- 
mogenization theory for such discrete elliptic operators (see for instance [7], [6]) ensures 
that there exist homogeneous and deterministic coefficients Ah om such that the solution 
operator of the continuum differential operator —V • ^4homV describes P-almost surely 
the large scale behavior of the solution operator of the discrete differential operator 
— V* • A(-,to)V. As for the periodic case, the definition of Ah om involves the so-called 
correctors <p : Z d x — > R, which are solutions (in a sense made precise below) to the 
equations 

- V* • A(x,uj)(^ + V(t>(x,uj)) = 0, xeZ d , (3.2) 

for £ G R d . The following lemma gives the existence and uniqueness of the corrector <fi. 
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Lemma 1 (corrector). Let f2, P, {Oz} z ^z<i, and A be as in Definitions® ® anc/[7| Then, 
for all £ G M. d , there exists a unique measurable function (ft : Z, d x $7 — > IR such that 
0(0, •) = 0, V(ft is stationary, (V</>) = 0, and (ft solves (|3.2p P-almost surely. Moreover, 
the symmetric homogenized matrix A^ om is characterized by 

Z ■ ^hom£ = ((£ + V0) • A(£ + V0)) . (3.3) 

As mentioned in the introduction, the standard proof of Lemma Q] makes use of the regu- 
larization of (|3.2|) by a zero-order term \x > 0: 



IMft^x,^) - V* • A(x,lj)(£ + V<ftJx,v)) 



0, 



x 6 



(3.4) 



Lemma 2 (modified corrector). Lei P, {#z} 26 z d > an d A be as in Definitions® ® anc/[?| 
Then, for all fi > an<i £ £ K , i/iere exists a unique stationary function 4>^ E L 2 (S7) which 
solves p.4p T-almost surely. 



In order to proceed as in the periodic case and use a spectral approach, one needs to 
suitably define an elliptic operator on L 2 ($7) (which is the stochastic counterpart to the 
space Hp CT (Q) of Section [2]). Stationarity is crucial here. Following [9], we introduce 
differential operators on L 2 (f2): for all u G L 2 (0), we set 



Bu( 



u(9 ei uj) — u{lo) 



u{0 ed )-u{oj) 



u(lo) — u(9- ei uj) 



U{U)) - u(6^ ed Uj) 



(3.5) 



We are in position to define the stochastic counterpart to the operator of Definition [2j 



Definition 8. Let Q, 

C : L 2 (Q) -> L 2 {tt) by 



{6> 2 } 2gZ d, and A be as in Definitions EJ El and [3 We define 



Cu(u) = -D*-i4(w)Du(w) 

= y^^ 0i2 (n(q;) - n(6> 2 a;)) 

2~0 



where D and D* are as in (13.51). 



In probabilistic terms, the operator — C is the generator of the Markov process called the 
"environment viewed by the particle". This process is defined to be (9x t w )> where (Xt) 
is a random walk whose jump rate from a; to a neighbor y is given by oj x ,y 

Using Definition [8] and the stationarity of (ft^, Lemma [2] implies that (ft^ is the unique 
solution in L 2 (Q) to the equation 

Ou + £)0 M = 0, (3.6) 

where 

:= D* -A(u)C (3.7) 

At the level of the corrector (ft itself (which is not stationary), the weak form of (13. 6f) 
survives for \x = 0: for every ^ S L 2 (fi), we have 

(D^- AD (ft) = (Dift-A£). (3.8) 
For all / £ L 2 (fi), we let £(/, /) be the Dirichlet form associated with C, defined by 

(Df.ADf) = ±J2 ( w o,*(/(0, w) - /H) 2 } . (3.9) 
16 
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As in the periodic case, the homogenized diffusion matrix satisfies the identity 

£-4hom£ = (£- AO -£(<f>,<}>). (3.10) 

The proof is formally the same as for (|2.2p . provided we use the weak form (|3.8p of the 
corrector equation, which holds for <p in place of ip (although 4> is not stationary). 

We refer the reader to [7] for the proofs of the statements above. 

3.2. Spectral representation and approximations of the homogenized coeffi- 
cients. The operator C is bounded, positive, and self-adjoint on L 2 (0). By the spectral 
theorem, for any function / £ L 2 (J7), we can define the spectral measure e/ of C projected 
on /, that is such that for any bounded continuous function ^ : R + — > 1R, one has 

</■*(£)/) = f *(A)de,(A). 

JM.+ 

As in the periodic case, we can express the homogenized diffusion matrix in terms of the 
spectral measure projected on Z). 

Lemma 3. Let f2, P, {0 z } Z £Z d > A, and C be as in Definitions [5l [?| anc/0. We let Ah om 
denote the associated homogenized diffusion matrix (j3.3j) . and d be the local drift (|3.7p . 
Then, the following identity holds 

S-A bam £={£-AQ- [ jde a (A), 
Jr+ a 

where is the spectral measure of C projected on d . 
Proof. In view of formula (|3.10p . we need to show that 

£(0,0) = I \ de a (A). 
Jr+ A 

This is either a consequence of Kipnis and Varadhan's arguments (see in particular [HI 
Theorem 8.1]), or a consequence of [H Corollary 1 &: Remark 2]. We detail the second 
argument. [H Corollary 1 & Remark 2] imply that lim M _j.o = V0 strongly in L 2 (Q), 
hence 

lim£(0 A1 ,0 M ) = £{4>,4>). 
Besides, for all fi > 0, we have by definition of the spectral decomposition 

and the result follows by the monotone convergence theorem. □ 

From Lemma O we deduce that the approximations A^ ^ introduced in Theorem [1] and 
further characterized in Proposition [2] may also be used in this discrete stochastic case, 
provided the notation (•) is understood as the expectation (instead of periodic average). 
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Err 



D,q,k 



(M) 



3.3. Suboptimal estimate of the systematic error. We let fl^fc, Pk, and A^k be as 

in Section [5J In order to quantify the systematic error, we introduce, for any D,q,k > 0, 
the function Erro^^ : R + — > R defined by 

/i 2fc i£k<D/4, 
/i 2fc ln+ +9 (/i- 1 ) ifJfe = D/4, 
H D / 2 \j§.{lT x ) i£k>D/A, 

where we write ln + (x) = maxjlnx, 1}. The purpose of this section is to show the following 
theorem. 

Theorem 3. Let ft, P, {0 z } Z £Z d > A, and C be as in Definitions^ [6j [?| and\^ and e a be 
as in Lemma We let Ah orn denote the associated homogenized diffusion matrix (j3.3[) , 
and A^ be the approximation (|2.22p of A-^ oxn for fj, > 0, and k > 1. Then, there exists 
q > (depending on a and (3) such that for all £ G M. d with |£| = 1, 

Err 2 , g ,fc(//) ifd = 2, 

Err dAfc (^) if5>d>2, 

Err 6) i ifc (/i) ifd = 6, 

Err 6i0 ,fc(/u) z/ 12 > d > 6, 

Err d _ 6i0 ,fc(//) if d> 12. 



< e • (A^fc - A hom )£ < 



In order to prove Theorem [31 we need to introduce some vocabulary. For all 7 > 1 and 
q > 0, we say that the spectral exponents of a function / G L 2 (fi) are at least (7, —5) if 
we have 

Jo 

Note that, if (7', —q') < (7,-9) for the lexicographical order, and if the spectral exponents 
of / are at least (7, —q), then they are at least (7', —<?')• Hence, the phrasing is consistent. 

In order to prove Theorem [31 we first express the systematic error in terms of the spectral 
exponents of d. This is the object of Theorem |U We then prove estimates on these 
exponents in Theorem [5j which concludes the proof of Theorem [3) 

Theorem 4. Within the notation and assumptions of Theorem^ the following two state- 
ments hold: for all £ G R. d with |£| = 1, 

(1) If the spectral exponents o/U are at least (7, —q), then 

if j > 2k + 1, 



< e • (\,k - A hom )t < 

(2) Conversely, 



I'- 
ll 



2 k 



In^ifi- 1 ) ifj = 2k + l, 
/u 7 " 1 InU/x- 1 ) i/7 < 2ib + l. 



£ • 04 M|jt - ^ h om)£ > M + z"" 1 / de (A). 

JO 

This theorem extends [H Proposition 9.1]. We begin by proving the following result. 



Lemma 4. If the spectral exponents ofd are at least (7,-5), i/ten 



< £ • (A 



4om)C < /X 2 * + H 



7-1 



2A- 



: ln 9 ((|iu) -1 ) du. 



(3.11) 
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Proof of Lemma^ First, recall that 

{ • - A„„)f = A(M+A)2 ..Vv + A)3 <MA). 

The integral of the r. h. s. is non-negative and bounded by 

1 



de 8 (A). (3.12) 

/k+ + A ) 

We perform a sort of integration by parts on this integral. To this aim, we let /'(A) be 
given by 

fm d 1 = ( / x + A) 2fc - 1 (^ + (2fc + l)A) 
J 1 j <9AA(/x + A) 2fc A 2 (// + A) 4fc 

We then rewrite the integral (|3.12p in terms of /', and use Fubini's theorem: 

1 



oo r+oo 

\=0 Js=x 

f'(S) I" de (A) dS. 

5=0 JX=0 

We split this double integral in two parts, and treat the cases 5 € (1, +oo) and 5 £ (0, 1] 
separately. We begin with the case when 5 ranges in (1, +oo). We bound the inner integral 

r5 roo 

/ de B (A) < / de (A) = (d 2 ) < 4/3 2 < 1, 
Jx=o Jx=o 

by definition of the projection of the spectral measure on D. This yields for the first part 
of the double integral 

r+oo t-5 i 

We now turn to the case when 5 ranges in (0, 1]. The assumption on the spectral exponents 
of c) implies 

rl rS rl 

f'(5) de (A) d5 < / f'(5)S^]n q (6~ 1 ) d5. (3.13) 

6=0 J\=0 JO 

Noting that 

f'(5) < (2k + 1) 



8 2 { l i + 5) 2k ' 

we bound the r. h. s. of f|3. 13|) by (2k + 1) times 

-1 



ln 9 ^" 1 ) d5. 



lo ^ + S) 2k 

A change of variables yields the announced result. □ 

Proof of part |I]) of Theorem^ We first assume that 7 > 2k + 1. In that case, we let 
7' be such that 2k + 1 < 7' < 7. Since the spectral exponents of D are at least (7',0), 
Lemma [4] ensures that 

< e • (a^ - A hom )a < f k + m 7 ' -1 [" ( f~\ 2k < M 2fe . 

JO K 1 + u ) 
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We now turn to the case when 7 < 2k + 1. We need to estimate the integral of the r. h. s. 
of (|3.1ip . To this aim, we note that 

ln^Oiu)- 1 ) = (ln^" 1 ) -ln(u)) 9 < 2 9 ( ln^/i" 1 ) + \ln{u)\ q ) , 

so that the integral in (|3.1ip may be estimated by 



,.7-2 



^-^(rn^) + |ln(^r)du< 



ln«+V -1 ) if 7 = 2fc + 1, 



ln 9 (iT 



if 7< 2Ar + l, 



as desired. 

Proof of part (0|) 0/ Theorem [^} Let 5 > be such that 

/ de a (A) > 0. 
J 

By the non-negativity of the spectrum and of the integrand, 



£ ' [Au.k ~ ^hom)^ 



> 



A(^ + A) 2 ---(2 fc ~V + A) 2 

2k(k-l)^2k rS 



5(fj, + 5) 2 • • • (2 k ~ V + <5) 2 
Hence, 

£ • {A^ k - A hom )£ > fi 2k . 
In addition, there exists C > such that for all A G (0, /x], one has 

,2fe C 



de (A) 
deo(A). 



A* 



A(/i + A) 2 ---(2 fc -V + A) 2 



> 



Therefore, 



£ • - Aom)e > A^ 1 / de (A), 

Jo 

which concludes the proof of the theorem. 

It remains to estimate the spectral exponents of 0. 



□ 



□ 



Theorem 5. Within the notation and assumptions of Theorem there exists q > 
depending only on the ellipticity constants a and (3 such that the spectral exponents of d 
are at least 

' (2,-q) ifd = 2, 

(d/2 + 1,0) ifb>d>2, 
(4,-1) ifd = G, 

(4,0) if 12 >d> 6, 

(d/2-2,0) if d> 12. 

Remark 1. We conjecture that the spectral exponents of D are in fact (d/2 + 1,0) for 
d > 2. If true, this would imply that the systematic error is in fact bounded by Err^o,*;^) 
for any d > 2 and k. 

In order to prove Theorem [5J we will make use of the following result. 
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Lemma 5. Within the notation and assumptions of Theorem [3J there exists q > de- 
pending only on the ellipticity constants a and f3 such that 

M - 2 ln"( fl - 1 ) i/rf = 2, 



where ^2 is as in Definition^ 



if5>d>2, 
l n +(» ifd = 6, 

1 ifd>6, 



Lemma [5] is a consequence of the results of [3] . Its proof, which is slightly technical, is 
deferred to Appendix lAl 

Proof of Theorem [31 For all A < /i, one has 

t >1 

( M + A) 2 (2/x + A) 2 " 36 

Hence, 

/ de (A) < 36/ / - — — —-2 de B (A). 
Jo Jr+ iP + A) 2 (2^ + A) 2 

The announced bounds then follow from Lemma [5] for d < 12. 

For d > 13, we use instead [H Theorems 2.3 and 2.4], which ensure that there exist C > 
such that for all /i > 0, 

Ta^cMA) < C/ /2 " 3 . 
jo 

This shows that the spectral exponents of d are at least (d/2 — 2,0), since 

/" M <MA) < /i X^de^X). 
Jo Jo 

□ 

3.4. Complete error analysis. As for the periodic case, <^ can be accurately replaced by 
4>n,Ri the solution of the modified corrector equation on a finite box Qr with homogeneous 
Dirichlet boundary conditions. We refer the reader to [2] for details. 

In order to perform a complete error estimate, one still needs to estimate the variance term 
in the r. h. s. of the identity corresponding to (ll.3|) . This is the object of the following 
theorem. 

Theorem 6. Let f2 ; ¥, {0 z } Z £z d > an< ^ A be as in Definitions^ [6], cmc?[?| We let j4hom 
denote the associated homogenized diffusion matrix (|3.3p . and for all k > 1, fj, > 0, and 
L > 0, we define the approximation Ank,L o/^hom as 

e-^w* : = «(^+v0 M )-A(e+v0 M ))>i 
fe-i it-ife-i 

+^ J] VkA^^L + l*J2Y^ U k,i,j((<t>2i»<i>2^))L, 
i=0 i=0 j>i 

where the coefficients % ^ and Vk,i,j are as i n Proposition the modified correctors 
are as in Lemma\^ and denotes the spatial average 

h i-> ((h)) L ■= / h(x)xL(x)dx, 
Jz d 
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where x i— > xl{x) is an averaging function on (—L,L) d such that f^ d XL(x)dx = 1 and 
W^Xl\\l°° ^ L~ d ~ x . Then, there exists an exponent q > depending only on a, f3 such 
that 

(L~ 2 + /x 2 )ln« /i" 1 if d = 2, 
L~ d + /i 2 L" d+2 i/ d > 2. 



Theorem [6] is a direct consequence of [31 Theorem 2.1 &; Remark 2.1] applied to each term 
of Auk i n the form f|2.25|) of Proposition [2j 

3.5. Polynomial decay of the variance along the semi-group. We end this section 
with a short remark concerning some results of [8]. Let (St)t>o be the semi-group associ- 
ated with the infinitesimal generator — C introduced in Definition [HJ In [Sj, the asymptotic 
decay to of the variance of Stf is investigated. A slight modification of [HJ Theorem 2.4] 
reads as follows. 

Theorem 7. Let f G L 2 {Vt) be such that (/) = 0, and let 7 > 1, q > 0. The following 
two statements are equivalent : 

(1) The spectral exponents of f are at least (7, —q) ; 
(2) 

<(^/) 2 ><^ 7 

From Theorem O we thus obtain the following result, which strengthens [HI Theorem 2.3 
and Corollary 9.3] when 4 < d < 12. 

Corollary 1. Within the notation and assumptions of Theorem 0. there exists q > 
depending only on the ellipticity constants a and (3 such that 

t~ 2 ln q + (t) ifd = 2, 
r (d/2+i) if 5 >d>2, 



((St*) 



< 



t" 4 ln+(i) ifd = 6, 



r 



t -(d/2-2) 



if 12 > d > 6, 
»/ d > 12. 
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Appendix A. Proof of Lemma [5] 



We adopt the notation of @J. In particular, we set T = fj, , denote by Gt the Green's 
function associated with the elliptic operator T~ 1 — V* • AV, 4>t the associated modified 
corrector, and we set tpT '■= Note that Gt and 4>t depend on the diffusion coefficients 
A. The claim of the lemma is equivalent to 

T 2 In 9 T 

rp3-d/2 

InT 
1 
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< 



if 
if 
if 
if 



d = 2, 
5>d>2, 
d = 6, 
d>6. 



(A.l) 
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Since {tpr) = 0, it holds that {ipj-) = var [tpr]- From the identity tJjt = T(cj)T — ^ 2 t) we 
learn that ipT depends continuously on the diffusion coefficients by [31 Lemma 2.6] so that 
one may apply the variance estimate of [31 Lemma 2.3]. In particular, 



var [ijj> 



T 



< 



sup 



dco P , 



(A.2) 



where the sum runs over the edges of Z d 
We proceed in four steps. 
Step 1. Proof of 



sup 

OJe 



duj P _ 



< 



(\VM*)\+»d(T)(l + \V<far(z)\))G T (0,e) 
+(l + |Vfer(z)|) I G T (0,w)G T (e,w)dw, 



(A.3) 



where e = (z, z + e$), Gt(0, e) := &r(0, 2 + e$) — Gt(0, z), Gt{z, w) = Gj-(z + ej, w) — 

Estimate (|A,3P is a direct consequence of [H 



1 if d > 2. 



G T (z,w), and // d (T) 
(3.10) & (3.21)], and [3 (2.14) & (2.16)]. 
S'tep 2. Proof of 

<(|WtW| 2 + Hd{T)\l + \V4>2t(z)\ 2 )G t (0, ef) dx 



< 



T 2 \n q T if d = 2, 

T if d = 3, 

InT if d = 4, 

1 if d>4, 



(A.4) 



where g only depends on the ellipticity constants a, (3. To prove (jA.4j) . we first replace 
the gradient of the Green's function by the Green's function itself and appeal to the 
deterministic optimal pointwise estimate of [31 Lemma 4]: 

|G T (x,e)| < G T (x,z) + G T (x,z + e t ) < fi d (T)(l + \x - z\f~ d min{l, VT \x - zf 1 }. 

By stationarity, (|VVt(^)| 2 ) = (|W>t(0)| 2 ), and {\V ct> 2T {z)\ 2 ) < 4 <|^ 2T (0)| 2 >, so that by 
[31 Proposition 2.1] and @j (3.27) & (3.29)], 

<(|VM*)| 2 + »d(T?(X + |V^ 2T (z)| 2 )G T (0, e) 2 > dx 

< «|V^rW| 2 > + / u d (T) 2 <|^ 2T (0)| 2 )) / / u d (T) 2 (l + |z|) 2 ( 2 " d )min{l,yT 2 |z|- 2 }dz 

li d (T)9(T+l)T if d = 2, 

(VT + 1)VT if d = 3, 

lnT+1 if d = 4, 

1 if d>4, 



< 



which yields (|A.4j) for T > 1. 



23 



24 



A. GLORIA & J.-C. MOURRAT 



< 



Step 3. Proof of 

(l + |V</> 2T (z)| 2 ) / / G T (0,w)G T (0,w')\G T (e,w)\\G T (e,w')\dwdw'^dz 

T 2 \n q T if d = 2, 
T 3-d/2 if 5 > d > 2, , 
InT if d = 6, [A ' b) 
1 if d > 6. 

where g only depends on the ellipticity constants a, /3. 

We first estimate the Green's function using the deterministic pointwise estimate of [31 
Lemma 4]: 

[ ((l + \V<p 2T (z)\ 2 ) [ [ G T (0,w)G T (0,w')\G T (e,w)\\G T (e,w')\dwdw')dz 
Jz d \ Jz d Jz d I 

< I I M^fa + kD^a + klf-Vinl^x/Tl^l^^minll,^^'!- 1 }* 1 
Jz d Jz d 

x / ((l + \V(b2T(z)\ 2 )\GT(e,w)\\G T (e,w')\)dzdwdw' (A.6) 
Jz d 

for some k > 1 to be fixed later (A; = 5 will be enough). We then deal with the inner 
integral, and appeal to the Meyers' estimate of [31 Lemma 2.9] and the bounds of [31 
Proposition 2.1] on the moments of the modified correctors. We let p > 2 be the Meyers' 
exponent. By Holder's inequality in probability with exponents ((p — 2)/p,2/p), Cauchy- 
Schwarz' inequality, and stationarity of VGt, we have 



< 



<(1 + \V<P2T(z)\ 2 )\G T (e,w)\\GT(e,w')\)dz 

(l + |V0 2T (z)| 2 ^- 2 ))(|G T ( e , U ;)n 1 / p <|G T (e, U ; / )l P > 1/Pd ^ 

< J (l + |V0 2T (z)| 2p/(p - 2) ) (\V z G T (z - w,0)\ p ) 1/p (\V z G T (z - w',0)\v) 1/p dz 

< n d (Ty [ (\V z G T (z-w,0)\ p ) 1/p (\V z G T (z-w',0)\ p ) 1/p dz. (A.7) 

Jz d 

The combination of (jA.6|) & (]A.7p with [31 Lemma 2.9] yields 

/ /(1 + |V0 2T (*)| 2 ) f [ G T (0,w)G T (0,w')\G T (e,w)\\G T (e,w')\dwdw')dz 
Jz d \ Jz d Jz d I 

< Hd(T) q / / gT(\w\)g T (\w'\)h T (z - w)h T (z - w')dzdwdw', 
Jz d Jz d Jz d 

where gr(t) = (1 + i) 2 ~ d min{l, VTt~ l } k , and /it is such that: for R ~ 1, 

/ M*) 2 < 1, 

J\x\<R 

and for all R 3> 1 and all j > 1, 



/ /i T (x) 2 dx < (2^' J R) d - 2 ( d - 1 )min{l, VT2^'} fe 

i2J_R<|x|<2J+ 1 J R 
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As we shall prove in the next step, this implies (|A.6h . Combined with Step 1 and Step 2, 
this proves the lemma. 



Step 4. Proof of 



gT(\w\)gT(\w \)Ht(z — w)Ht{z — w')dzdwdw' 



< 



T 2 InT if d = 2, 

T 3-d/2 if 5 > > 2, 

InT if d = 6, 

1 if d > 6. 



(A.f 



The proof of (|A.8|) is made technical because the bounds on Ht do hold integrated on 
dyadic annuli, and not pointwise. In line with the bounds on hx, we prove the claim by 
using a doubly dyadic decomposition of Z d x Z d combined with the results of [31 Proof 
of Lemma 2.10, Steps 1, 2 & 4], that we recall for the reader's convenience: there exists 
R ~ 1 such that for alii G N, 



/hx{z)hT{z — x)dzdx 
t z\<\z-x\ 

(2 J i?) 2 max{l,ln(v / T(2 J J R)" 1 )} if d = 2, 
(2 i R) 2 if d>2, 



(A.9) 



x\<4R J \z\<\z— x\ 



tiT{z)hT(z — x)dzdx 



InT if d = 2, 
1 if d > 2. 



(A.10) 



We first use the symmetry with respect to w and w' to restrict the set of integration to 
\w'\ > | w |, and we make a change of variables 



9T{\w\)gT{\w'\)hT(z — w)hx{z — w')dzdw' dw 



< 2 / g T (\w\)g T (\w'\)h T (z - w)h T (z - w')dzdw'dw 

Jwez d Jw'ez d ,\w'\>\w\ Jz d 

= 2 I f [ g T (\w\)g T (\w'\)h T (z)h T {z-(w-w'))dzdw'dw, 

JweZ d J w-w'£Z d , \w'\>\w\ Jz d 
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followed by the associated doubly dyadic decomposition of space 
<?t(|H)5t(|w'|)/it(z — w)hr(z — w')dzdw'dw 
~ /oiR^i 'i <-9.+i B 9t{\w\)9t{\w'\) 

x / h T (z)h T (z - (w - w'))dzdw'dw (A.ll) 

^ /i ,^ti 2 3 R<\w-w' <2 3+1 R ^ Vl l/y Vl ly 
jeN J\w\<RJ [w , UH 

x / h T (z)h T (z -(w- w'))dzdw'dw (A. 12) 



x/ h T (z)h T (z- (w -w'))dzdw'dw (A. 13) 



uj|<R 



+ / /,. „. D «r(|to|)fl!r(Kl) 



' - w'\ < R 
v | > \w\ 



x / h T (z)h T {z - (w - w'))dzdw'dw, (A.14) 



where ~ 1 is as above. We begin with the last term (|A.14|) of the sum, and appeal to 
(jA.lOp and the definition of qt'- 



w\<R 



w -w'\<r 9t{\w\)9t{\w'\) / h T (z)h T (z - (w - w'))dzdw'dw 

w'\ > \w\ 



< / / h T {z)h T {z - {w - w'))dzdw' 

J \w-w'\<R . 



< 



lnT if d = 2, . . 

1 if d>2. 



We continue with ()A,13p , Since \w — w'\ < R, gxdw'l) ~ <?t(M), and we have using 
(jA.lOp and the definition of gx'- 

V/ / 9t{\w\)9t{\w'\) \ h T (z)h T (z - (w - w'))dzdw'dw 

< ( V / g T (\w\) 2 dw) I [ [ h T (z)h T (z - (w - w'))dzdw') 

\~7Z J2 i R<\w\<2 i + 1 R / \ J\w-w'\<R JlA I 



< 



'2 i fl<|w|<2 i + 1 .R J \J\w-w'\<R 

TlnT if d = 2, 

VT if d = 3, ( . s 

lnT if d = 4, (A - i5j 
1 if d>4. 
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For (1X12]) we note that \w\ < R and 2 j R < \w-w'\ < 2 j+1 R imply that \w'\ ~ 2 j R, and 
we appeal to (|A.9p : 



< 



< 



*\r»\<RJ ??<\ w ,- w '\<* +1 * 
g T (0)R d y2g T (^R) 

j£N J2JR<\w-w'\<Z'+ 1 R 

TlnT if d = 2, 

VT if d = 3, 

InT if d = 4, 

1 if d > 4. 



Pt(|w|)pt(^'|) / h T (z)h T (z - (w - w'))dzdw'dw 



hT(z)tiT(z — (w — w'))dzdw' 



(A.17) 



The dominant term is (|A~TT|) . We split the double sum into three parts according to the 
range of z and j: 

• the diagonal part: \i — j\ < 1, 

• the off-diagonal parts: i> j + 2 and j > i + 2. 

For |z — j| < 1, we use the inequality \w\ < \w'\ < + \w — w'\ so that for the i th term 
of the sum, \w'\ ~ 2 l R. In particular, using (IA.9j) . this yields for the diagonal term 



E E 

*6N | i _ i |<l-'2 , *<M<2' 



< 



< 



^(2 i J R)V(2 i J R) 2 



x / /i T (z)/i T (2; — (u; — -u/))d2idy/dw; 

hT(z)tiT(z — (w — w'))dzdw' 



2 i - 1 R<\w-w'\<2 i + 2 R 



]T(2*i?) 



d+2(2-d) 



min{l,VT(2 J J R)" i } & (2 i J R)V d (T) 



< 



/i,(T)^(2 i J R) 6 - d min{l,VT(2' J R)- 1 } E 

igN 

T 2 InT if d = 2, 

VT 3 if d = 3, 

T if d = 4, 

\/T if d = 5, 

InT if d = 6, 

1 if d > 6. 



(A.18) 



We turn to the first off-diagonal term: those integers i,j such that i > j + 2. In this case, 
we use the estimate \w — w'\ — \ w\ < \w'\ < \w — w'\ -\-\w\, which shows that for the (i,j) th 
term of the sum, \w'\ ~ 2 l R. In particular, using (lA.9j) . this yields for the first off-diagonal 
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term 



EE f f 



5T(|w|)g T (k'|) 



^— ' ^— ' L D .i 1^0,4.1 D / 2 J il < lie - tu'l < 2 i+1 il 

x / hT{z)hx{z — (w — w'))dzdw'dw 



< 



iGN j<«-2 



2i_R<|«;-u;'|<2i+ 1 i? 



hT(z)tiT(z — (w — w'))dzdw' 



< ^(2^)^(2-,) min{1 ^ |2 ^ r i }5 J- (2iR) 2 pi d (T) 

< Y^{^R) d+2{2 ~ d) min{l, \/T|2 i J Rr 1 } 5 (2 i i?)Vd(T) 

igN 

T 2 lnT if d = 2, 

^/T 3 if d = 3, 

T if d = 4, 

if d = 5, 

lnT if d = 6, 

1 if d > 6. 



(A.19) 



We now treat the last term of the sum, that is those integers i,j such that j > i + 2. 
Then, similarly to (lA.19j) we deduce that for (i,j) th term of the sum, \w'\ ~ 2 J R. Hence, 
using (IA.9D . we obtain 

/ 

J 2 3 R < \w - w'\ < 2 J + 1 _R 



E E 

imj > i+2 J*R<\ W \<2^RJ pj- 



< 



i>j+2 



E(2' i J R)V(2 4 J R) E 9T&R) j 
E ( 9T^R) [ 



g T (\w\)g T {\w'\) 

x / hT(z)h,T(z — (w — w'))dzdw'dw 

hr(z)hr{z — (w — w'))dzdw' 



2JR<\w-w'\<2i+ 1 R 



l2iR<\w~w'\<2J+ 1 R. 

< Y,9T{2 ] R){^Rf^ d (T) E WR) d g T (?R) 

J'SN «<J-2 

T 2 lnT if d = 2, 

^3 



I h T {z)h T {z -{w- w'))dzdw' I V {^R) d 9T(2 l R) 



i<7'-2 



< 



VT if d = 3, 

T if d = 4, 

VT if d = 5, 

lnT if d = 6, 

1 if d > 6. 



(A.20) 



as for the first off-diagonal term. 

Estimate fA8|) then follows from the combination of (fA3D - (fA~T4|) with ([ATT8l) - ([A~20l) 
and d3jSI)-(|3j2|). 
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Figure 1. Periodic cell in the discrete case 



Appendix B. Numerical tests in the discrete periodic case 

Numerical tests of [T] have confirmed the sharpness of Theorem [2] for the approximation 
A)j,,i,r,l on a discrete periodic example. In the present work, we consider the same discrete 
example, and numerically check the asymptotic convergence of A^^,r,l to ^4hom- As 
expected, the systematic error is reduced, and the limiting factor rapidly becomes the 
machine precision. The discrete corrector equation we consider is 

-V* ■ A^ + Vcp) = inZ 2 , 

where V and V* are as in (I3.ip . and 

A(x) := diag [uj XtX+ei , uj x , x+e2 ] . 

The matrix A is [0, 4) 2 -periodic, and sketched on a periodic cell on Figured! In the example 
considered, ui x ^ x + ei and ui x ^ x + e2 represent the conductivities 1 or 100 of the horizontal edge 
(x, x+e±) and the vertical edge (x, x + &2) respectively, according to the colors on Figure Q3 
The homogenization theory for such discrete elliptic operators is similar to the continuous 
case (see for instance |10j in the two-dimensional case dealt with here). By symmetry 
arguments, the homogenized matrix associated with A is a multiple of the identity. It can 
be evaluated numerically (note that we do not make any other error than the machine 
precision). Its numerical value is ^Ihom = 26.240099009901 .... To illustrate Theorem [2] 
in its discrete version (which is similar, see [T] for related arguments), we have conducted 
a series of tests for A^- I n particular, we have taken /i ~ i?- 3 / 2 , L = R/3, and a filter of 
infinite order. In this case, the convergence rate is expected to be of order 3 for A^i, and 
of order 6 for A^z. This is indeed the case, as can be seen on Figure El where R denotes 
the number of periodic cells and ranges from 6 to 400 (that is log(i?) up to 2.6). 
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